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ABSTRACT 


A.n  algorithm  is  described  for  generating  two-dimensional,  visible- 
line  projections  of  three-dimensional  objects  that  are  bounded  by 
patches  of  quadric  surfaces. 

The  main  task  of  the  algorithm  is  the  calculation  of  intersections 
between  quadric  surfaces.  A parameterization  scheme  is  used.  Each 
quadric-surface  intersection  curve  (QSIC)  is  represented  as  a set  of 
coefficients  and  parameter  limits.  Each  value  of  the  parameter  repre- 
sents at  most  two  points,  and  these  may  easily  be  distinguished.  This 
scheme  can  find  the  coordinates  of  points  of  even  quartic  (fourth-order) 
intersection  curves,  using  equations  of  no  more  than  second  order. 

Methods  of  parameterization  for  each  type  of  QSIC  are  discussed, 
as  well  as  the  problems  of  surface  bounding  and  hidden-surface  removal. 
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I . INTRODUCTION 


1.1  General 

This  report  describes  an  algorithm  for  the  computer  generation  of 
orthographic  as  well  as  perspective  projections  of  three-dimensional 
objects  whose  surfaces  are  made  up  of  segments  of  quadric  surfaces,  also 
known  as  quadric  patches. 

There  are  now  many  algorithms  for  generating  pictures  of  solid 
objects  whose  surfaces  consist  of  planar  polygons.  A polygon  is  a seg- 
ment of  a first-order  surface  (plane)  bounded  by  other  first-order 
surfaces . 

The  Braid  algorithm  [Braid  1975]  extends  this  to  allow  also  cylin- 
ders and  parts  of  cylinders.-  Cylinders  are  quadric  surfaces. 

This  paper  is  based  on  the  Woon  algorithm  [Woon  1970  and  Wocn- 
Freeman  1971] , which  processes  second-order  (quadric)  surfaces  bounded 
by  second-order  surfaces.  Quadric  surfaces  include  such  objects  as 
spheres,  cylinders,  cones,  ellipsoids  and  hyperbolic  paraboloids.  They 
are  not  too  complex  mathematically. 

An  important  advantage  of  using  quadric  patches  is  that  many 
objects  (especially  man-made  ones)  can  be  more  precisely  modelled  by  a 
small  number  of  quadric  patches  than  by  a large  set  of  small  polygons. 

Another  advantage  is  that  shaded  pictures  can  be  made  with  better 
control  over  such  effects  as  mach  banding  and  specular  reflection 
[Phong  1975]. 
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Previous  quadric  patch  algorithms,  such  as  [Mahl  1972]  and 
[Weiss  1966] , as  veil  as  the  Woon  algorithm,  become  laborious  when 
solving  for  the  intersection  of  two  quadric  surfaces.  This  is  a fourth- 
order  problem,  and  there  is  no  easy  method  for  finding  the  roots  of  a 
fourth-order  equation,  or  even  merely  determining  whether  there  are  any 
real  roots. 

The  algorithm  presented  here  reduces  each  fourth-order  problem  into 
a second-order  problem.  The  algorithm  uses  a parametric  representation. 
Every  point  on  the  intersection  of  two  quadric  surfaces  is  determined 
uniquely  by  a numerical  parameter  and,  in  many  (but  not  all)  cases,  by 
a Boolean  sign  parameter.  Since  the  fourth-order  quadric  surface  inter- 
section curves  (QSICs)  can  be  represented  parametrically,  the  computer 
memory  need  only  store  about  a score  of  numbers  to  represent  a QSIC, 
rather  than  storing  some  lengthy  point -by-point  description  for  each  QSIC. 

1 . 2 The  Woon  Algorithm 

Since  the  Woon  algorithm  is  not  well  known,  and  since  it  forms  a 
basis  of  this  algorithm,  the  Woon  algorithm  will  be  presented  first. 

?or  each  object,  the  input  data  for  the  Woon  algorithm  consists  of 
the  surface  equations,  surface  bounds,  surface  intersections,  and 
vertices.  Once  an  object  is  specified,  the  user  may  ask  for  different 
views  of  the  object,  representing  different  locations  of  the  vantage 
point.  Pictures  may  be  generated  in  either  orthographic  or  perspective 
projection,  and  with  the  hidden  lines  displayed,  dashed,  or  suppressed. 
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For  each  surface  the  user  must  specify  the  following: 

a)  The  equation  of  the  surface.  This  is  expressed  in  terns  of  the  ten 
coefficients  of  the  equation: 

q(x,y,z)  = q x2  + q y2  + q z2  + q xy  + q yz 

1 2 3 4 5 

+ qzx  + qx  + qy  + az  + q =0  (l) 

*6  7 8 ‘9  0 

The  quadric  surface  divides  the  three-dimensional  space  into  two 
open  regions:  the  exterior,  in  which  q(x,y,z)>0;  and  the  interior , in 

which  q(x,y,z)<0.  On  the  surface  itself,  q(x,y,z}=0. 

b)  The  bounds  of  the  surface.  When  using  the  Woon  algorithm,  one  is 
normally  interested  in  only  one  connected  segment  of  each  surface.  This 
surface  segment  (the  "quadric  patch")  is  defined  as  the  locus  of  points 
on  the  surface  satisfying  a certain  Boolean  condition,  which  is  specified 
for  each  surface. 

Each  bound  consists  of  another  quadric  surface,  and  a polarity 
( + or  -).  If  the  polarity  is  positive  (+),  the  bound  is  satisfied  (true) 
only  for  points  in  the  exterior  of  the  bounding  surface.  If  the  polarity 
is  negative  (-),  the  bound  is  satisfied  only  for  points  in  the  interior 
of  the  bounding  surface. 

Each  set  of  bounds  consists  of  one  or  more  bounds.  The  set  of 
bounds  is  satisfied  if  and  only  if  ail  bounds  in  the  set  are  satisfied. 
Each  patch  has  one  or  more  sets  of  bounds.  The  patch  is  the  locus  of  all 
points  on  the  surface  satisfying  at  least  one  set  of  bounds. 

c ) The  intersection  of  the  surface  patch  with  other  surface  catches. 

Woon  distinguishes  between  planar  surface  interesections  which  lie  in  a 
plane,  and  non-planar  intersections  which  do  not.  He  also  distinguishes 
between  protrusive  and  recessive  intersections. 
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In  addition  to  the  above  information,  needed  for  each  surface  and 
each  surface  intersection,  the  user  must  also  specify  all  real  vertices. 
These  are  points  where  three  or  more  intersection  curves  meet. 

Quadric  surfaces  have  the  property,  when  viewed  from  certain  view- 
points, of  "folding  in  back  of  themselves",  such  that  part  of  the 
quadric  surface  hides  another  part.  The  curve  at  which  this  happens, 
the  locus  of  points  of  tangency  of  rays  from  the  vantage  point  to  the 
surface,  is  called  a "limb".  (Astronomers  use  this  term  for  the  visual 
edge  of  a celestial  body. ) 

The  drawing  produced  by  the  Woon  algorithm  consists  of  edges. 

There  are  two  types:  real  edges,  which  are  portions  of  the  intersections 
of  quadric  surfaces;  virtual  edges,  which  are  portions  of  limbs. 

Woon  used  a variation  of  the  Loutrel  algorithm  [Loutrel  1970]  to 
determine  whether  or  not  points  are  visible.  This  requires  that  the 
quadric  patches  be  broken  up  into  sections  called  faces . A patch  which 
folds  in  back  of  itself  is  divided  into  two  or  more  faces  by  the  limb. 

A face  is  classified  as  a front  face  or  a back  face,  depending 
upon  its  orientation  with  respect  to  a viewer.  Only  front  faces  can  be 
visible.  The  orientation  of  the  faces  forming  an  edge,  as  well  as  the 
protrusiveness  or  recessiveness  of  the  edge,  help  determine  whether  the 
edge  is  potentially  visible. 

The  computation  of  pictures  is  divided  into  two  parts.  First,  for 
each  object,  the  object-dependent  quantities  are  computed.  These 
include  real  edges,  real  vertices,  and  the  partial  classification  of 


edges . 


Then,  for  each  view  of  the  object,  the  virtual  edges  and  virtual 
vertices  (intersections  of  virtual  edges  and  real  edges)  are  computed, 
and  finally  the  visible  portions  of  each  edge  are  found.  To  generate 
the  data  for  drawing  the  edges,  the  object  must  be  rotated  about  three 
angles:  azimuth  (<}>),  elevation  (0),  and  twist  (f).  For  perspective 

drawings,  the  distance  from  the  origin  to  the  viewpoint  (D)  and  the 
distance  from  the  viewpoint  to  the  picture  plane  (d)  must  be  considered. 

Woon  suggests  two  methods  for  calculating  points  on  a quadric  sur- 
face intersection  curve  (QSIC).  If  the  QSIC  is  planar  it  will  be  a 
conic  section,  and  will  be  relatively  easy  to  handle  mathematically. 

If  the  QSIC  is  not  planar,  then  the  curve  may  be  computed  as  a series 
of  points,  each  within  a specified  distance  from  its  neighbors.  The 
term  "vector"  is  used  for  the  directed  line  segment  from  one  of  these 
points  to  the  next. 

Each  vector  has  a dominant  direction,  either  ±x,  ±y,  or  ±z.  For 
the  first  vector  in  a curve,  a guess  is  made  as  to  its  dominant  direc- 
tion, and,  if  it  takes  the  curve  out-of-bounds  for  its  surfaces,  then  it 
is  rejected  and  another  direction  is  selected.  For  subsequent  vectors, 
a strong  hint  as  to  the  dominant  direction  is  provided  by  the  direction 
of  the  last  vector.  Newton's  method  is  then  used  to  obtain  the  best 
approximation  to  the  actual  direction. 

Woon  discusses  a method  of  determining  whether  an  intersection  is 
planar  or  not.  This  method  is  discussed  and  expanded  upon  later  in  this 
paper. 

The  Loutrel  algorithm  [Loutrel  1970]  is  used  by  Woon  to  determine 
whether  or  not  an  edge  or  point  is  visible.  Each  point  on  an  edge  has 
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an  order  of  invisibility,  which  is  the  number  of  front-faces  between  it 
and  the  vantage  point.  The  order  of  invisibility  changes  whenever  the 
line  from  the  point  to  the  vantage  point  intersects  an  edge  of  the  hid- 
ing surface. 

Loutrel's  original  algorithm  dealt  only  with  polygonal  faces.  All 
edges  were  straight-line  segments;  computing  the  intersecting  points  of 
edge  projections  was  exceedingly  simple.  This  is  not  the  case  for  the 
projections  of  edges  of  quadric  patches.  These  edges  may  be  as  high  as 
fourth-order  curves.  Woon  specifies  a very  complicated  algorithm  for 
determining  these  intersection  points,  as  well  as  an  algorithm  for 
propagating  the  order  of  invisibility  from  edge  to  edge.  Hone  of  this 
will  be  used  in  this  report. 

1.3  QUADRA'*' 

Woon  wrote  a program  called  QUADRAW  which  implements  a variation  of 
his  algorithm.  For  each  object,  QUADRAW  requires  information  on  scaling, 
vector  length,  and  tolerances. 

For  each  surface,  QUADRAW  requires  the  following: 

(1)  The  equation  of  the  surface.  (There  are  also  some  "auxiliary  sur- 
faces" which  are  not  part  of  the  object,  but  bound  other  surfaces.) 

(2)  The  bounds  of  the  surface. 

For  each  surface  intersection,  it  requires: 

(1)  The  two  intersecting  surfaces. 

(2)  The  type  (recessive  or  obtrusive). 

(3)  The  planarity  of  the  intersecting  surfaces  and  of  the  surface 


intersection. 


{U)  If  neither  surface  is  planar,  but  the  intersection  is,  the  equa- 
tion of  the  plane. 

(5)  Specification  of  each  real  vertex. 

QUADRAW  computes  all  curves  in  a piece-wise  linear  fashion.  Each 
QSIC  is  represented  in  memory  by  a series  of  small  vectors.  This 
approach  requires  use  of  a computer  with  a large  primary  memory. 

QUADRAW  does  not  use  Woon's  variation  of  the  Loutrel  algorithm  to 
its  fullest  extent.  It  does  use  classification  of  edges  to  determine 
which  edges  might  be  visible , but  then  uses  a brute-force  method  to 
compute  edge  visibility.  This  brute-force  method  works  nicely,  and  is 
used  here,  with  some  modification. 

As  each  view  is  processed,  the  maxima  and  minima  in  the  view-plane 
are  computed  for  each  surface.  For  each  potentially  visible  edge, 
visibility  is  computed  one  point  at  a time.  Each  point  is  checked 
against  each  face  whose  extrema  include  the  coordinates  of  the  point. 

If  at  least  one  face  is  found  which  hides  the  point,  then  the  point  is 
invisible . 

The  testing  to  see  whether  a face  hides  a point  on  the  edge  (test 
point ) proceeds  as  follows : 

The  line  of  sight  is  computed  from  the  vantage  point  to  the  tesu 
point.  The  equation  of  the  surface  is  reduced  by  considering  only  points 
along  the  line  of  sight.  This  results  in  a second-order  equation  is  most 
cases,  with  first-order  equations  sometimes  occurring.  If  the  solution 
of  the  equation  yields  two  roots,  one  corresponds  to  a point  on  a back 
face  and  may  be  discarded.  If  the  remaining  point  is  within  bounds,  then 
it  hides  the  test  point , and  no  further  test  is  needed  for  that  point . 
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II.  THE  QUADRATIC  FORM 


The  quadratic  surface,  as  represented  by  eq.  (l)  is  a quadratic 
fora  [Woon  1970,  Dresden  1930,  Newman-Sproull  1973].  It  may  be  repre- 
sented in  vector-matrix  form,  as  follows: 

Suppose  the  location  of  a point  in  3-space  is  represented  by  the 
vector:  x = (x  y z l).  Eq.  (l)  then  may  be  represented  as  follows: 

q(x,y,z)  = q(x)  = (x  y z l)i 
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or , for  short , 


(3) 


q(x)  = x Q x' 

where  Q is  the  U x U matrix  in  equation  (2).  Q is  called  the 
discriminant  of  the  quadric  surface.  Notice  that  Q is  symmetric  and 
that,  for  any  real  non-zero  scalar  8,  8Q  is  equivalent  to  Q,  since  they 
describe  the  same  surface.  In  this  report  the  same  symbol  will  on 
occasion  be  used  for  both  a quadric  surface  as  well  as  its  discriminant. 
An  alternate  form  which  will  also  be  used  is  the  following: 


A D F G 
D B E H 
F E C J 
_ G H J K 

such  that 

q(x,y,z ) = Ax2  + By2+  Cz2  + 2Dxy  + 2Eyz  + 2 Fzx 


(M 


+ 2Gx  + 2Hy  + 2Jz  + X = 0 


(5) 


■■ I 
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The  upper-left  3x3  principal  sub-matrix  of  the  discriminant  con- 
tains all  the  second-order  terms.  This  will  be  called  the  "sub- 
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The  rank  of  the  discriminant  and  the  sub-discriminant  are  helpful 
in  classifying  quadric  surfaces.  These  are  invariant  under  the  trans- 
formation described  below.  The  rank  of  the  discriminant  is  invariant  under 
any  nonsingular  transformation,  and  that  of  the  sub-discriminant  is  invariant 
under  any  transformation  in  which  the  upper-left  3x3  sub-matrix  is 
nonsingular. 


Transformations  of  the  Discriminant 


In  setting  up  the  description  of  an  object,  the  untransformed  space 
will  be  the  u-v-w  space,  and  the  transformed  space  will  be  the  x-y-z 
space.  The  transformation  matrices  will  be  F and  T,  each  being  the 
inverse  of  the  other. 

First  let  us  consider  the  trar.fcrmaticn  x = uF,  where  x = (x  y z 1) 
and  u = (u  v w 1),  1st 
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where  F is  a congruence  transformation.  The  sub-discrininant , F , is 


an  orthogonal  transformation.  If  is  orthonormal,  then  the  determin- 
ants of  both  F and  F are  unity, 
u 

Assume  that  “F  = F \ T is  the  transpose  of  F if  they  are  ortho- 

u u 

normal.  The  lower-right  element  of  T (the  constant  term)  is  unity. 


Suppose  the  equation  for  the  quadric  surface  P is: 

->  -Vp 

u P u =0  in  untransformed  space,  and 
x x =0  in  transformed  space. 

Using  the  transformation  x = uF  on  equation  (8b),  we  have: 

— ► — ^P 

u F "?  F u . Comparing  this  with  equation  (8a),  we  have 
T 

P = F ? F 


(3a) 

(8b) 


3y  pre-  and  post-multiplication,  we  have: 

T?r 


?F??*rT  = I?I=?,  or  T = 7 ? *FT 


(9) 


(10) 


It  is  possible  to  break  down  the  transformation  into  two  parts: 
rotation  and  translation.  They  may  be  applied  repeatedly  in  any  order. 
The  rotation  matrix  is: 
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The  sub-dsicriminant  is  affected  only  by  rotation.  Rotation  does 
not  affect  the  constant  term,  in  the  lower-right  corner  of  the  discrim- 
inant . 
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The  translation  matrix  is: 

T = 

Translation  does  not  affect  the  sub-discriminant.  The  constant 
term  is  affected  only  by  translation. 

A transformation  can  be  formed  by  a rotation  followed  by  a trans- 
lation (F  = RT),  or  by  a translation  followed  by  a rotation  (F  = TR) . 

T and  R do  not  usually  commute,  just  as  rotation  and  translation  to  not 
usually  commute. 

2.2  Canonical  Form 

It  is  possible  to  put  the  discriminant  of  a surface  into  a canoni- 
cal form,  in  which  the  axes  of  the  coordinate  system  are  the  axes  of  the 
object,  with  non-zero  diagonal  elements  of  the  discriminant  being  given 
precedence.  In  finding  the  canonical  form,  one  also  finds  the  transform- 
ation matrix  which  will  transform  the  discriminant  into  its  canonical 
form. 

In  the  following  algorithm,  the  discriminant  P is  transformed  into 
its  canonical  form  Q by  transformation  matrix  7.  The  matrix  E is  used 
as  an  intermediate  transformation.  The  operation  means  assignment, 
as  in  a computer  program. 

a)  First,  a rotation  is  performed  to  diagonalize  P . This  rotation 

rp 

is  represented  as  the  orthonormal  matrix  7 . We  have  Q *■  7 P F~ , 

u u u u u 
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If  has  both  non-zero  and  zero  eigenvalues,  the  non-zero  ones  should 
be  placed  towards  the  upper-left  corner  of  Q , and  the  zero  ones  towards 
the  lower-right . 

*?  is  then  forced  by  augmenting 


«-  '0 

ru  ! o 

'o 

0 0 0 1 


We  have  Q, 


and  Q is  of  the  form 


A 0 0 G 
0 3 0 H 
0 0 C J 
G H J K 


b)  We  must  now  try  to  eliminate  as  many  cff-diagonal  elements  as  possi- 
ble. (Skip  this  step  if  A = 3 = C = 0.)  We  take  Z as: 


10  0 0 


0 0 10 

6 6 5 1 

_x  y 


6 = -G/A; 

x 

6^  = -H/3  if  3 i 0, 

5 = -T/C  i?  C r 0, 
z 


6 = 0 if  3 = 0; 

y 

5 = 0 if  C = 0. 


We  next  have  the  "standard"  transformation: 


P Z P Z\ 


■>-  l T. 


(15a, b) 


■ 1 | 1 — ■ — -™-— - m 


r 


- ik  - 

c)  At  this  point  P is  in  one  of  four  forms 
If  P is  in  the  fom: 

[' 

then  it  is  already  in  canonical  form 
If  it  is  in  the  form: 

i 

~A  0 0 o" 

0 3 0 0 
0 0 0 J 

I 

00  JK 


A 0 0 0 
0 B 0 0 
0 0 C 0 
0 0 0 X 


then  it  is  a paraboloid.  See  (d),  below,  on  how  to  eliminate  the  K 
(constant)  term. 

If  it  is  in  one  of  these  forms: 


'aooo* 

~0  0 0 G~ 

0 0 0 K 

0 0 0 H 

0 0 0 J 

0 0 0 J 

_0  H J X _ 

_G  H J K_ 

then  we  must  remove  all  off-diagonal  pairs  except  for  the  top-  (left-) 
most  one.  This  may  be  done  as  follows: 

For  the  case  on  the  left,  one  uses  the  transformation  matrix: 


" 1 

0 

0 

-i 

0 

0 

K 

s 

0 

(16) 

0 

-3 

K 

0 

_ o 

0 

0 

]_ 

k 
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with  K2  + S2  = I 


, C H 

and  — = — 

o J 


This  removes  one  of  the  off-diagonal  pairs,  making  J=0.  One  then 
uses  the  standard  transformation  (I5a,b). 

For  the  case  on  the  right,  this  must  be  done  twice,  once  to  remove 
H,  and  once  to  remove  J. 

d)  At  this  point,  P has  one  off-diagonal  element.  If  K = 0,  we  have 
the  canonical  form.  If  not,  we  make  a translation  of  the  same  form  as 


eq.  (lU)  except  with  the  following  6’s: 


o = -K/2G  if  G 4 0, 
x 

6 = -K/2H  if  H 4 0, 

y 

5 = -2C/J  if  J / 0, 

z 


6 = 0 if  a = 0; 
x 

6 = 0 if  H = 0; 

y 

6 = 0 if  J = 0. 


After  applying  the  standard  transformation  (I5a,b),  we  have  the  final 
canonical  form  for  P,  and  the  final  transformation  matrix  7. 


2.3  The  Discriminant  Form  for  Conic  Sections 

The  conic  sections  are  quadric  curves,  being  also  of  the  quadratic 
form.  The  discriminants  of  a conic  are  3x3  matrices,  and  the  sub- 
discriminants are  the  2x2  upper-left  subnatrices  of  the  discriminants. 
Other  than  the  reduction  in  the  number  of  dimensions,  the  conic  sections 
are  completely  analogous  to  quadric  surfaces;  the  transformation  matrices 
are  slightly  simpler: 


f f 0 

ux  uy 

f f 0 

vx  vy 

6 6 1 

x y 


K S 0 

-3  K 0 

5 6 1 

x y 


o o 

where  X“  + S 
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2 . U Classification  of  Quadric  Surfaces  and  Conic  Sections 

Tables  I and  II  are  guides  to  the  classification  of  quadric  sur- 
faces and  conic  sections,  respectively.  They  are  adapted,  in  part,  from 
a chart  given  by  Dresden  [Dresden  1930,  p.  230]. 

One  of  the  columns  in  these  tables  is  the  absolute  value  of  a 
signature . The  signature  of  a matrix  is  the  number  of  positive  eigen- 
values minus  the  number  of  negative  eigenvalues. 

In  Table  I,  we  use  the  following: 


*11  q12  ql3  qll 


-21  q22  -23  q2U 

q31  q32  q33  q3l+ 

ql*l  a-l;2  qU3  qhU 


A D F G 
D 3 E H 
F E C J 
G H J K 


A D F i 


D 3 E 
F E C 


d = Rank(Q); 


m = Rank ( Q ) ; 


s = abs( Signature (Q_  )), 
u 


3 

l a.. 

tii  *1X 


= A + 3 + C 


3 3 <1, , q . , 

l l ° 

1-1  J-i+1  qJi  4jj 


= A3  + AC  + 3C  - D2  - Ez  - F2  (22) 


3 h q.  . q 
y y ii  ij 

i=l  J=i+1  qji  qjj 


T + T K - G2  - H2  - J2 


f 
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SINGULAR  SURFACES 


d 

m 

s 

Conditions 

Peal  surface  (Imaginary  mart  in  parens.) 

Planar 

X 

0 

0 

INVALID 

1 

X 

Coincident  Planes 

2 

0 

0 

Single  Plane 

2 

1 

1 

D2>0 

INVALID  (Imag.  Parallel  Planes) 

2 

1 

1 

d2<o 

Two  Parallel  Planes 

2 

2 

0 

t2<o 

Two  Intersecting  Planes 

2 

2 

2 

t2>o 

LINE  (Two  Intersecting  Imaginary  Planes) 

Non-planar  Singular 


3 

1 

1 

Parabolic  Cylinder 

3 

2 

0 

t2<o 

Hyperbolic  Cylinder 

3 

2 

2 

T2>0;  TjD3<0 

Elliptic  Cylinder 

3 

2 

2 

T2>0;  TjD3>0 

INVALID  (Imag.  Cylinder) 

3 

3 

1 

a 

Cone 

3 

3 

3 

3 

POINT  (Imaginary  Cone) 

NON-SI 

NGU 

LAP. 

SURFACES  d=L ; 

m=2  if  det (Q  ) = 0 
u 

0 der.o 

tes 

zer 

•o; 

m=3  if  det(Q^)  ? 0 

+ denotes  positive; 


denotes  negative; 


± denotes  nonzero. 


(Sj  IGI  S 

0+0 
0-2 
± + 1 

± - 1 

± + 3 

± - 3 


Conditions 

(T2<0) 

(t2>o) 

a 

a 

3 

3 


Surface 

Hyperbolic  Paraboloid 
Elliptic  Paraboloid 
Hyperboloid  of  One  Sheet 
Hyperboloid  of  Two  Sheets 
INVALID  (Imaginary  Ellipsoid) 
Ellipsoid 


Conditions:  a: 
6: 


T >0;  det(Q  ) x T <0;  or:  T <0. 
z u * 

T2>0;  det(Qu)  x Tj>0. 


Table  I.  Guide  for  the  Classification  of  Quadric  Surplus 
(adapted  from  Dresden  1930) 


- 19 


r*u 

qi2 

%3  _ 

~A 

D 

n 

\x  \ 

r\  _ 

o'  — 

ro 

C-22 

q23 

= 

D 

3 

H 

; Qu  = 

r, 

H, 

— 

D 

G 

K 

K 

D 

E 

n 

L 31 

C\J 

m 

C7* 

q33  _ 

— 

— 

t 

d = Rank(Q);  a = Rank(Q  ); 

u 


= a, , + a„  = A + 3 
l *xl  ~et2 


s = ab 3 ( Signature ( Q )) 
u 


D 


2 3 

l l 


i=l  j=i+l 


= A3  + 3X  + AX  - D2  - G2  - H2 


ji  a s_  Conditions 

Singular 

10  0 
ill 
2 0 0 

2 1 1 D > 0 

2 

2 1 1 D,  < 0 

2 2 0 IQ  I < 0 

U 

2 2 2 j Q | > 0 


Curve 


Iir/ALID 

Coincident  Lines 
Single  Line 

INVALID  (Imag.  Parallel  Lines) 
Two  Parallel  Lines 
Two  Intersecting  Lines 
POINT  (Inter.  Imag.  Lines) 


3 

1 

1 

Parabola 

3 

2 

0 

ISJ 

< 

0 

Hyperbola 

3 

2 

2 

ISJ 

> 

0; 

|Q|Tj  < 0 

Ellipse 

3 

2 

2 

> 

0; 

JO 

V 

o 

INVALID  (Imaginary  Ellipse) 

Table  I-.  Guide  for  the  Classification  of  Ccnic  Sections 
(adapted  from  Dresden  1930) 


r 
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III.  CLASSIFICATION  OF  QUADRIC  SURFACE  INTERSECTIONS 

Given  two  quadric  surfaces,  the  first  two  questions  that  arise  are: 

1)  Do  they  intersect?  and 

2)  If  so,  what  is  the  nature  of  their  quadric  surface  intersection 
curve  (QSIC )? 

3 . 1 The  Pencil  of  Two  Quadric  Surfaces 

Suppose  we  have  two  quadric  surfaces,  with  discriminants  ? and  Q. 

—y  — ►rri 

In  matrix  form,  the  equations  for  the  two  surfaces  are  x P x*  = 0 

-e  -tf? 

and  x Q x*  = 0. 

The  equation 

x (Q-a?)xr  =0  (26) 

represents,  for  all  real  values  of  a( finite  or  infinite),  a surface  on 
the  "pencil"  of  ? and  Q. 

For  a = 0,  we  have  the  surface  Q.  For  a = ± 00 , we  have  the  surface 


If  surfaces  P and  Q intersect,  then  their  intersection  (QSIC)  is 
the  "’case  curve"  of  the  pencil,  and  it  lies  in  all  the  surfaces  of  the 
pencil.  This  "base  curve"  is  not  to  he  confused  with  the  "base  curve" 
of  a parameterization  surface,  which  will  be  used  extensively  later  on 
If  the  two  surfaces  do  not  intersect,  then  none  of  the  real  sur- 
faces of  the  pencil  intersect.  In  addition,  the  pencil  contains  some 
imaginary  surfaces,  and  these  may  be  among  those  listed  as  INVALID  in 
Table  I. 


* Since  SP  represents  the  same  surface  as  ? for  any  real  r.on-sero  scalar 
3;  if  3=1 /ot,  lim  (1/a) (Q-aP)  = lim  (Q/a)  - ? = -F,  which  is  equiva- 


The  general  form  fcr  a member  of  the  pencil  of  P and  Q is  given  by 


S(a)  = Q - aP.  (2 

3.2  Classification  of  Pencils 

One  may  classify  pencils  according  to  the  classification  of  the 
simplest  surface  in  the  pencil.  If  the  pencil  has  a base  curve,  then 
the  base  curve  (a  QSIC)  has  the  same  classification  as  its  pencil. 

If,  for  some  a,  R( a)  has  rank  one  or  two,  then  P.(a)  represents 
either  a plane  or  pair  of  planes.  This  pencil  is  called  "planar". 

Failing  this,  if,  for  some  a,  R ( Cc ) has  rank  three , the  pencil  is 
"non-planar  singular". 

If  P.(a)  is  never  singular,  then  it  always  has  rank  four , and  the 
pencil  is  non-singular . 

Two  surfaces  do  not  intersect  if  one  of  the  following  occurs: 

1)  For  some  value  of  ct,  R(a)  is  INVALID;  or 

2)  Zither  P or  Q does  net  intersect  some  R(a). 

3.3  Determining  Whether  a Surface  Intersection  is  Planar  or  Singular 
If,  for  some  value  of  a,  R(a)  has  rank  of  two  or  less,  then  the 

pencil  is  planar.  As  Woon  has  pointed  out  [woon  1970,  pp.  3^-36],  this 
occurs  when  two  conditions  are  met  for  the  same  value  of  ct: 

(1)  Det(R(a))  = 0;  ana 

(2)  The  sum  of  the  3x3  principal  minors  of  R(ct)  vanishes.  This 

may  be  expressed  as:  D3(P.(ct))  = 0. 

If  condition  (l)  is  met  but  condition  (2)  is  not,  then  the  surface 
is  non-planar  singular. 
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IV.  THEORY  OF  PARAMETERIZATION 

In  this  section,  methods  of  parameterization  are  discussed.  The 
parameterization  is  done  in  a u-v-w  coordinate  system.  A congruence 
transformation  is  then  used  to  transform  point  coordinates  into  the 
x-y-z  system.  The  parameter  is  denoted  by  t. 

1.1  Parabola 

Suppose  we  have  a parabola  of  the  form  Au2  + 2Hv  = 0.  Taking 
m = -A/2H,  we  have  v = mu2.  The  parameter  equations  are 

u = t ; v = mt2 . (28) 


1* . 2 Ellipse 

Suppose  we  have  an  ellipse  of  form  Au2  + Bv2  + K = 0.  Taking 

r = /-K/A  and  r = v-'fj 3 as  being  the  semi-axes,  we  have 
u v 


u * v 

The  parametric  equations  are 


u 


and 


(29) 


This  form  is  to  be  used  only  -l<t<+l,  in  which  case  v will  have 

only  positive  values.  (A  more  complete  method  is  given  in  section  k.k.) 

This  parameterization  is  well-behaved  for  values  of  t within  the 

range  [-1,  +i].  That  is,  if  ds  = /duz  + avz  , ds/dt  does  not  vary  too 

much  if  r and  r are  of  the  same  order  of  magnitude, 
u v 


- 2h 


b.3  Hyperbola 

Suppose  we  have  a hyperbola  of  the  form  Au2  + Bv2  + K = 0. 


Taking  r = /K/A  and 
u 


r = /-K/B,  ve  have 
v 


u = 


73^  V 


l+tz 

V = 1 IT  r ! 

1-t  v 


(30) 


for  -l<t<+l. 


A hyperbola  has  two  disjoint  parts.  This  parameterization  gives 
only  one  of  them  the  one  for  which  v>0. 

In  contrast  to  the  parameterization  of  the  ellipse,  this  form  is 
not  well-behaved.  At  values  of  t approaching  ±1,  small  changes  in  t 
result  in  large  changes  in  u and  v.  However,  the  hyperbola  is  then  very 
close  to  its  asymptotes  and  is  practically  a straight  line. 

4,4  Mors  ConcLe^e  Ferns 

~:r  cer.tr-..  s :•  Licses  and  hyperbolae ) , the  following  modifi- 
cations may  be  u3ec: 


Select  a parameter  t’  which  takes  :n  the  values  from  -2  and  +2 
inclusive  for  ellipse--.  / r - /r e - : iae,  too  values  -2,  0,  and  +2  are 
excluded. ) 


For  t’  > 0, 
For  t*  < 0, 
Thi3  yields  for 

u * 

and  hyperbolae 

u = 


we  have  t = t ' 
we  have  t = t ' 
illicses 


o* 


- 1 and  a = +1 ; 
+ 1 and  a = -1 . 


(31) 


(32) 
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Since  for  an  ellipse, u and  v should  he  periodic  functions,  we  can 
set  up  a parameter  t",  which  is  normalized  to  t'  hy  adding  or  subtract- 
ing multiples  of  four,  such  that  t'  is  in  the  range  [-2,  +2]. 

U. 5 Non-nlanar  Intersections 

The  parameterization  of  planar  intersections  is  described  above. 

For  non-planar  intersections,  the  curve  lies  in  a quadric  surface  which 
has  a "base  curve"  which  is  either  a line,  a parabola,  a hyperbola,  or  an 
ellipse.  This  base  curve  siay  be  a cross-section  of  the  quadric  surface, 
in  a plane  perpendicular  to  the  main  axis.  There  is  one  set  of 
straight  lines,  (one  line  passing  thru  each  point  of  the  base  curve) 
such  that  each  line  lies  wholly  in  the  quadric  surface,  ana  every  point 
on  the  quadric  surface  lies  on  one  of  these  lines. 

One  uses  parameterization  to  select  a point  on  the  base  curve  and 
its  corresponding  line.  3y  solving  a quadratic  (second-order)  equation, 
one  can  then  find  the  intersection  of  the  line  with  any  other  quadric 
surface.  In  this  manner,  all  the  points  of  the  QSIC(s)  may  be  found. 

SPECIAL  NOTE : THE  FOLLOWING  ALGORITHM  HAS  NOT  BEEN  IMPLEMENTED , AND  IS 

PRESENTED  ONLY  AS  A GUIDE  TO  A FUTURE  IMPLEMENTATION . 


V.  DATA  NEEDED  3Y  THE  ALGORITHM 


The  algorithm  needs  the  following  input  information  in  some  form: 

(1)  For  each  object,  it  needs  to  know  the  dimensions  of  the  "object 
box",  which  is  a cube  or  rectangular  parailepiped  in  which  the  object 
is  contained,  as  well  as  the  resolution  and  vector  lengths  which  will 
be  used. 

(2)  Surface  ecuatior.s  must  be  known  for  each  surface.  These  may  be 
presented  either  as  ten  coefficients  (in  the  form  of  eq.  (l)),  or 
presented  as  follows: 

a)  The  type  of  surface  (ellipsoid,  hyperboloid,  cylinder,  etc.). 

b)  The  lengths  of  radii,  semi-axes,  etc. 

c)  The  orientation  and  displacement  from  the  origin. 

An  interactive  graphics  terminal  could  be  used  for  specifying  a 
surface  and  then  manipulating  and  distorting  it  to  suit. 

(3)  Rounds  must  be  specified  for  each  quadric  patch  in  the  same  manr.e 
as  in  the  Wocn  algorithm. 

(L)  Surface  intersections  must  be  specified.  For  each  intersection 
the  user  must  specify  to  the  algorithm  the  two  intersecting  patches, 
the  multiplicity  of  the  intersection  (how  many  disjoint  parts  it  has), 
and  whether  it  is  to  be  a "smooth"  or  a "sharp"  intersection. 

A "sharp"  intersection  abruptly  separates  two  surfaces,  providing 
a clear  edge  which  can  be  seen  if  viewed  from  a proper  angle. 

"Smooth"  intersections  are  used  when  several  quadric  patches  are 
used  to  approximate  a single,  higher-order  surface.  These  are  not 


— 


m 
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included  in  drawings  unless  they  occur  along  limbs.  Often  the  first 
derivatives  will  be  continuous  across  a smooth  intersection.  If  one 
wished  to  approximate  a torus  (donut)  by  using  patches  of  ellipsoids, 
hyperboloids  of  one  sheet,  and  cones,  smooth  intersections  would  be 
used. 

One  usually  specifies  intersecting  surfaces  as  bounds  for  each 
other.  However,  if  there  is  a smooth  intersection  with  a continuous 
first  derivative,  then  it  may  be  difficult  to  tell  whether  a point  is 
on  one  side  of  a boundary  or  the  other.  Therefore,  the  program  should 
automatically  compute,  for  every  smooth  intersection,  another  surface 
from  the  pencil  of  the  two  intersecting  surfaces.  This  surface  will 
act  as  an  auxiliary  bounding  surface, 
ing  surfaces  at  a sharp  angle. 


It  should  meet  both  intersec 
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VI.  CLASSIFYING  SURFACE  INTERSECTIONS 


Consider  two  intersecting  surfaces  P and  Q.  If  either  is  planar, 
then  the  intersection  is  planar.  The  algorithm  described  in  section 
VII  will  apply  in  these  cases. 

If  neither  of  the  two  intersecting  surfaces  is  planar,  then  one 
must  find  the  "simplest"  surface  of  the  pencil,  according  to  the  list 
of  Table  III.  The  simplest  surfaces  are  those  in  which  the  subdiscrim- 
inant is  singular.  We  first  solve  the  equation: 


det(R  (a))  = det(Q  -aP  ) = 0.  (33! 

u u u 

In  order  to  insure  that  this  has  at  least  one  real  root,  we  will  dis- 
allow the  case  where  jP  | is  singular  and  |Q^J  is  non-singular.  (If 
this  should  happen,  we  interchange  them.)  Equation  (33)  is  of  (at  most) 
third  order.  This  may  be  written  as: 


-det(P  )ct3  + K ct2  - K.a  + aet(Q  ) = 0, 

U 2 * U. 


with  K being  the  sum  of  the  determinants  of  the  combinations  of  two 

columns  of  ? and  one  of  Q , and  E being  the  sum  of  the  determinants  of  the 
u u l 

combinations  of  one  column  of  ? and  two  columns  of  0 . 

u u 

Once  a single  root  is  found,  synthetic  division  may  be  used  to  find 


the  other  real  roots,  if  they  exist. 

As  pointed  out  in  [Woon  1970,  pp.  3^-36]  R(a)  is  planar  when 
det(R(a))  = 0 and  the  sum  of  the  3x3  principal  minors  of  R(a)  also 
vanishes.  (This  last  condition  may  be  expressed  as:  D,(R(ct))  = C.) 
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These  conditions  may  only  occur  when  aet(Ru(cO)  = 0.  Therefore,  each 
root  in  a is  tested  to  see  whether  R(a)  is  planar.  If  any  R(a)  are 
planar,  the  simplest  is  then  chosen  as  a parameterization  surface. 

If  none  of  the  R(a)  so  far  tested  is  planar,  one  checks  to  see 
whether  any  one  is  a hyperbolic  paraboloid.  If  one  is,  then  it  is  used 
for  parameterization  (see  Section  8.U).  If  not,  then  one  of  them  may 
be  a cylinder,  which  would  then  be  used  for  parameterization. 

If,  at  any  point,  an  INVALID  R(a)  is  found,  then  the  two  surfaces 
do  not  intersect. 

Table  III  contains  a list  of  the  surfaces  which  may  be  used  for 
parameterization.  The  simplest  surfaces  are  near  the  top  of  the  list. 


L 
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VII.  PLANAR  INTERSECTIONS 


A pencil  is  planar,  if,  for  some  a,  det(R(a))  = 0,  ana 


D, (R(a) ) = 0. 


7.1  Finding  and  Classifying  Intersections 

First  cr.e  must  determine  the  nature  of  the  intersection  (according 
to  Table  I)  and  transform  the  planar  surface  into  canonical  form.  Then 
one  selects  some  surface  of  the  pencil,  other  than  R(cc),  and  applies 
the  same  transofrmation  to  it.  This  surface  will  usually  be  either  ? 


a)  If  the  intersection  lies  in  a LINS,  then,  in  the  space  of  the  can- 
onical form,  this  line  is  the  w-axis.  In  this  case,  one  cancels  the 


first  two  rows  and  the  first  two  columns  of  the  discriminant 

j Q r ^ 

ether  surface,  getting  a 2 x 2 matrix  of  the  form  jtj  . 

If  C - J = K = 0,  then  the  line  is  the  intersection. 

If  C = J = 0 and  X r 0,  then  there  is  no  intersection. 
Otherwise,  taking  *D  = J2-CK,  we  have: 


No  intersection  if  J‘-CK  < 0; 


ijra  mterse 


ectior.  point  if  J2-CK  = 0; 


Two  intersection  points  if  J2-CX  > 0. 
b)  If  the  intersection  lies  in  a single  plane,  or  in  two  coincident 
planes,  the  plane  is  the  u-v  plane  in  the  space  of  the  canonical  form. 
One  takes  the  discriminant  of  the  other  surface,  cancels  out  the  third 
row  and  third  column,  and  compares  the  resultant  3x3  conic  section 
discriminant  against  Table  II  to  find  the  nature  of  the  intersection. 


_f  tr.ere  os  no  intersection. 


J 
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c)  If  the  intersection  lies  in  two  planes,  one  must  factor 
intersection  into  its  two  constituent  planes.  If  there  are 
planes,  we  have  an  equation  Au2  + K = 0,  with  AK  < 0,  which 
into  /j A j u ± /[iiy  = C. 


the  diplar.a 
two  paralle 
is  factored 


If  there  are  two  intersecting  planes,  we  have  Au2  + 3v2  = 0,  with 
A3  < 0,  which  is  factored  into  /] A f u ± /j'B]  v = 0. 

Then,  each  of  these  planes  must  be  put  into  single-plane  canonical 
form , making  the  appropriate  transformation  to  the  ether  surface  for 
each  plane.  Finally  each  plane  must  be  processed  separately,  as  above. 

If  it  should  happen  that  a planar  intersection  consists  of  two 
straight  lines,  these  must  be  processed  separately. 

7 . 2 Parameterization 


A f*  w ST* 

(unless  the 
One  then  rc 
conic.  For 
hyperbola, 
make  approp 
transformat 


the  nature  of  the  planar  intersection  is  found,  one  has 
intersection  was  a LINZ),  a conic  section  in  the  u-v-piane. 
tates  the  u-v  plane  so  that  the  axes  are  the  axes  of  the 
a parabola,  the  v-axis  is  the  axis  of  symmetry.  For  a 
the  u-axis  separates  the  two  parts.  One  must,  as  always, 
riate  transformation  of  both  the  other  surface  and  of  the 
ion  matrix,  *F.  One  then  applies  the  parameterization,  as 


outlined  ir.  section  17. 

When  actually  tracing  points,  one  is  in  x-y-z  space,  not  u-v-w 
space.  Therefore,  it  would  be  more  efficient  to  get  the  x-y-z 
coordinates  of  a point  directly,  without  having  to  transform  coordinate 
each  time,  and  without  taking  up  primary  memory  space  to  store  the 
transformation  matrix.  The  form  outlined  below  accomplishes  just  this. 
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X = 

(a  t 

X 

2 + b t 

X 

+ c ) / 6 + d 
X X 

y = 

(a  t 

y 

2 + b t 

y 

+ c ) /o  + d 

y y 

•7  — 

(a  "t, 
z 

2 + b t 
2 

+ c ) / 6 + d ^ 

where 

6 = 

i 

for 

a parabola  o: 

6 = 

i + 

2 

t : or 

an  ellipse; 

6 = 

1 - 

t2  for 

a hyperbcla. 

Using 

this 

form. 

one  may  perform  she 

setting  up 

the  c 

oeff  i 

cier.ts  ( 

a , b , c . d) . 

For  completeness,  one  nay  use  the  forms  outlined  in  section  U.U. 
Table  IV  contains  a FORTRAN  routine  to  do  this. 


*7  ”3  TT ■?  m ~ "O' 


rinding  rar 


o y»o*ri  a* 


er  Limits 


Not  all  of  a Quadric  Surface  Intersection  lurve  (Q3IC)  is  part  of 
an  edge.  Some  or  all  of  a QSIC  may  be  out  of  bounds  for  one  or  both 
of  the  intersecting  surfaces.  Since,  for  planar  intersections,  each 
value  of  the  parameter  represents  a unique  point,  the  endpoints  of  an 
edge  may  be  represented  as  parameter  values. 

laeh  QSIC  has  a multiplicity,  which  is  the  number  cf  distinct  edges 
which  are  intersection  curves  of  the  same  two  surface  patches.  An 
example  of  this  would  be  two  cones,  back-to-back,  with  a thin  circular 
cylinder  dividing  their  intersection  into  two  parts.  This  intersection, 
between  the  two  cones,  has  a multiplicity  of  two . 

The  intersection  cf  a sphere  with  a cylinder  of  diameter  less  than 
that  of  the  sphere,  such  that  the  center  of  the  sphere  is  along  the  axis 
cf  the  cylinder,  is  not  a single  intersection  of  multiplicity  two,  but 


_ 


J 


\n 
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Mote:  In  the  following  program,  x,  y,  and  z are  X(l),  X(2),  and  X ( 3 ) ; 

a^  is  A(l,x),  b^  is  A(2,x),  is  A(3,x),  and  d^  is  A(4,x); 
t is  T,  and  t'  and  t"  are  TP. 

SUBROUTINE  PLANAR  (TP,  A,  ITYPE,  X) 

C TP  IS  THE  PARAMETER,  A IS  THE  COEFF.  MATRIX,  h. 

C X IS  THE  RETURNED  CO-ORDINATE  LOCATION. 

C ITYPE  = -1  HYPERBOLA 

C 0 PARABOLA  OR  LINE 

C +1  ELLIPSE 

REAL  A(U,3)  X(3) 

IF  (ITYPE)  2,  3,  1 

ce:jtral  conics 

1 TP  = TP  - 1*#INT(TP+3IGN(2,TP)  )/4) 

2 T = TP  - SIGN(l . ,TP) 

DENOM  = SIGN  (l.  + ITYPE*-*T , TP) 

GO  TO  4 

C—  PARABOLA  OR  LINE 

3 T = TP 
DENOM  = 1. 

COMPUTE 

4 DO  5 I = 1,  3 

5 X(I)  = (A(1,I)»T#T+A(2,I)*T+A(3,I))/DEN0M  + A(4,l) 

RETURN 

END 


j 


! 


Table  IV.  FORTRAN  routine  for  solution  of  eq.  (35) 


- 35  - 


Figure  1:  one  for  each 

Intersection. 

rather  two  separate  intersections  of  multiplicity  one,  as  the  two  parts 
of  the  cylinder  are  two  separate  patches  (see  Fig.  l). 

A QSIC  will  have  one  pair  of  parameter  limits  for  each  unit  of 
multiplicity. 

If  a quadric  surface  intersection  is  a point  or  pair  of  points, 
then  any  points  which  are  out -of -hounds  must  be  rejected.  Otherwise, 
parameter  limits  must  be  found.  The  program  must  find  one  or  more  ranges 
of  the  parameter  for  which  the  parameter  is  continuously  within  bounds 
for  both  of  the  intersecting  surfaces. 

First  the  program  establishes  a list  of  bounding  surfaces  for  the 
QSIC.  Then  a scan  is  made  of  the  possible  values  of  the  parameter, 
searching  for  a value  whose  point  is  within  bounds. 

For  an  ellipse,  the  parameter  is  taken  from  -2  to  +2  (in  t'  form), 
with  parameter  increments  of  about  0.1.  For  a parabola  or  hyperbola, 
the  limits  are  the  values  of  the  parameter  at  the  limits  of  the  object 
box,  but  the  increments  for  the  hyperbola  should  be  about  0.1  also. 
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As  this  is  done,  the  value  of  q^(x)  = x for  each  bounding 

surface  ^ is  saved.  If  no  point  satisfying  the  bounds  is  found  in  the 

initial  scan,  then  one  may  look  for  "dips"*in  the  values  of  the 

a, (x)'s.  This  way  one  may  find  a small  intersection  which  was  over- 
b 

looked  before.  Interval  halving  is  performed  at  each  dip  to  see 
whether  it  results  in  a previously  missed  small  intersection. 

Once  a point  within  bounds  is  found,  one  can  begin  looking  for  the 
edge's  endpoints.  Basically,  this  is  tracing  the  edge,  for  different 
values  of  the  parameter,  until  a point  out  of  bounds  is  found.  (Often 
the  first  point  within  bounds  has  an  already-traced  neighboring  point 
out  of  bounds.)  When  such  a pair  of  points  is  found,  one  may  use 
interval  halving  to  locate  the  precise  end  point.  If  it  is  known  which 
bounding  surface  function  changes  sign  to  cause  one  point  to  be  in 
bounds  and  the  other  out  of  bounds,  one  may  use  Newton's  method  on  that 
"critical"  bounding  surface  function. 

If  one  is  looking  for  another  edge  on  a particular  plane  and  cannot 
find  one  within  the  parameter  values  not  included  in  the  range  of  the 
first  edge,  it  is  possible  that  the  first  edge  is  really  both  edges, 
with  a small  gap  in  between.  Then  one  should  go  back,  checking  all 
points  where  the  value  of  any  bounding  surface  function  changes  sign, 
to  see  whether  these  points  are  endpoints  of  such  a small  gap.  If 
none  is  found,  it  is  likely  that  the  user  made  an  error. 

Special  note  on  ellipses:  If  there  are  no  bounds  on  an  ellipse 

that  cause  parameter  limits,  then  the  limits  are  <-2,+2>.  If  the 
point  t'  = -2(+2)  is  in  the  middle  of  an  edge,  using  the  t"  format 
allows  one  to  extend  paramater  limits  above  +2  or  below  -2, 


An  extremum  towards  zero. 


A 


- 37  - 

VIII.  NON-PLAiJAR  QSICS 

A QSIC  not  contained  in  a plane  is  non-planar  or  twisted. 

Appendix  A shows  that  such  a QSIC  is  always  located  in  a quadric  surface 
with  at  least  one  set  of  non-intersecting  straight  lines  passing  through 
every  point  on  the  surface. 

This  section  shows  how  to  utilize  this  property  to  parameterize  a 
non-planar  QSIC.  Given  a parameter,  one  may  find  a point  (or  pair  of 
points)  by  solving  an  equation  of  no  more  than  second-degree.  Since 
the  quadratic  equations  is  easy  to  solve,  this  method  gives  the  exact 
location  of  each  point.  When  there  are  two  points  for  one  parameter 
value,  then  one  is  associated  with  an  increasing  parameter,  and  the 
other  with  a decreasing  parameter. 

8.1  Basic  Techniques 

A non-planar  QSIC  is  contained  in  a quadric  surface  which  shall  be 
called  the  "parameterization  surface".  It  is  one  of  the  surfaces  of 
the  pencil  of  the  two  intersecting  surfaces.  Table  III  contains  a 
list  of  these. 

Selecting  the  R(ot)  which  is  "easiest"  to  handle,  we  find  the  proper 
parameterization  for  its  base  curve.  For  a cylinder,  the  base  curve  is 
a cross-section  of  the  parameterization  surface,  in  a plane  perpendi- 
cular to  the  main  axis. 

Given  a parameter  one  finds  the  corresponding  point  on  the  base 

curve  using  eq's.  (32).  These  coordinates  are  referred  to  as  x , y , 

o o 

and  zfl . 


1 
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One  then  selects  any  member  of  the  pencil  other  than  the  para- 
meterization surface.  Usually  this  will  be  P or  Q.  One  finds  the 
straight  line  in  the  parameterization  surface  which  passes  through 
point  (xfl,  yQ , zfl),  and  solves  for  the  intersection  of  the  line  and  the 
other  surface.  This  will  be  a second-order  equation. 

Suppose  Y , Y > and  y are  the  direction  cosines  of  the  line 
x y z 

corresponding  to  the  base  curve  point  xQ  , yQ , zg.  The  location  of  any 
point  on  the  line  is  given  by 
x = x + Y s 

0 X 

y = y0  + Yys  (36) 

z = Z0  + Yzs 

where  s is  a secondary  parameter.  If  the  other  surface  is  surface  Q 
and  has  an  equation  of  the  form  of  eq.  (l),  then  the  second-order  equa- 
tion in  s is  of  the  form: 

a s2  + b s + c =0  (37) 


where 


a = qY2+qY2+qY2+<lYY  + qY  Y + q Y Y 
1 x 2 y 32  4 x'y  5 y 2 6 z x 

b = 2q  x Y + 2q  y y + 2q  z y + q (x  Y +y  Y„)  + q (y  Y,+2  y ) 
1 0 x 2 0 y 3 0 z 4 0 y 0 x 5 0 z 0 y 

+ q (z  Yv  + X y ) + q y + q y + q y 
6ox  oz  7X  8 y 9Z 


c=qx2+qy2+qz2+qxy  + qyz  + qzx  + qx 

10  20  30  400  500  600  70 


q z + q = q(x  ,y  ,z  ] 

9 0 0 0 0 0 


+ q y + 
a 0 
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The  quadratic  discriminant  (not  to  he  confused  with  the  discriminant 
matrix  of  a surface)  is  *D  = b2  - Uac.  If  *D  < 0,  there  are  no  real 
roots  in  s,  and  the  corresponding  value  of  t is  actually  invalid. 

If  *D  = 0,  then  there  is  one  root  in  s.  If  V > 0,  s has  two  real 
roots.  There  are  as  many  intersection  points  along  the  line  as  there 
are  roots  in  s.  When  tracing  a QSIC,  one  uses  the  larger  value  of  s 
when  the  parameter  is  increasing,  and  the  smaller  value  of  s when  the 
parameter  is  decreasing.  This  is  discussed  in  more  detail  in  section  IX. 

8.2  Cylindric  Intersections 

For  cylindric  intersections,  the  line  through  each  point  is  parallel 
to  the  axis  of  the  cylinder,  which  is  the  w-axis.  One  uses  the  substi- 
tution 


with  (e  e e ) s ( 0 0 1)  7 , 
x y z u 

The  cross  section  of  a cylindric  intersection  may  be  taken  anywhere 
along  the  main  axis  but  the  u-v  plane  is  simplest. 
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8.3  OSICs  Lying  in  Hyperbolic  Paraboloids 

A hyperbolic  paraboloid  has  two  sets  of  straight  lines  which  lie 
totally  in  the  surface.  These  are  called  reguli.  Each  regulus  is  a 
set  of  mutually  skew  lines,  and  each  point  on  the  paraboloid  is  the 
intersection  of  one  line  from  each  regulus. 

To  find  the  coefficients  of  parameterization,  one  takes  the  hyper- 
bolic paraboloid  P and  the  other  surface  Q,  and  transforms  them  both  by 
a congruence  transformation  so  that  P is  in  canonical  form: 

A 0 0 

0 -B  0 

P ~ 0 0 0 
0 0 J 


0 

0 


0 / 


(38) 


If  A is  negative,  reverse  the  sign  of  every  element  in  P. 

Now  transform  both  P and  Q by  the  following  valid  transformation, 
which  is  a product  of  a rotation  and  a scaling  transformation.  It  will 
distort  the  surfaces  but  not  change  their  nature: 


/T72A  -JTm  0 

/l72A  /l72B  0 


0 

0 


\ 


0 

0 


0 -1/J  0 

0 0 1 


The  result  is  that  Q is  in  the  form 


Q = 


- hi  - 


This  is  equivalent  of  the  equation  (in  u-v-w  space)  uv=w.  If  t is  the 
primary  parameter  and  s is  the  secondary  parameter,  u=t,  v=s,  and  w=st. 

In  the  x-y-z  coordinate  system,  we  may  express  the  coordinates  in 
terms  of  equations  (36)  as  such: 

x = (ct+d)  + (e  + f t)s 

XX  XX 

y = (ct+d)  + (e  +f  t)s  (39) 

J y y y y 

z a (ct+d)  + (e  + f t)s 

z z z z 

If  *F  is  the  transformation  matrix  from  u-v-w  space  to  x-y-z  space, 
then  the  co-efficients  are: 


(UO) 


t 
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IX.  PARAMETER  LIMITS  AND  QSIC  TRACING 


9.1  Basic  Parameter  Limits 

For  planar  QSICs,  parameter  limits  are  defined  in  a two-part  pro- 
cess. In  the  first  part,  the  basic  limits  are  found.  For  an  ellipse, 
the  basic  limit  string  is  <-2.  ,+2.>.  For  hyperbolae  and  parabolae, 
which  extend  to  infinity,  the  basic  parameter  limits  are  those  defined 
by  points  where  the  curve  meets  the  limit  of  the  object  box. 

Once  basic  limits  are  found,  the  bounded  parameter  limits  are 
found,  using  the  process  outlined  in  section  7-3. 

For  non-planar  intersections,  finding  the  basic  parameter  limits 
is  more  complicated.  These  limits  are  the  limits  of  the  range  of  values 
of  the  parameter  such  that  both  the  following  conditions  are  met: 

a)  If  the  base  curve  is  a line,  a parabola,  or  a hyperbola,  the  line 
corresponding  to  the  parameter  must  not  lie  totally  outside  of  the 
object  box;  and 

b)  Equation  (37)  has  at  least  one  real  root.  This  is  equivalent  to 
saying  that  the  quadratic  discriminant  V = b2-l+ac  is  non-negative.  One 
would  find  those  ranges  of  the  parameter  t'  for  which  V >_  0.  The  limits 
of  these  ranges  will  be  values  of  t'  for  V = 0.  Of  course,  if  for  all 
values  of  the  parameter  t'  satisfying  condition  (a),  V < 0,  then  there 
is  no  intersection. 

Once  the  basic  limits  are  found,  basic  limit  strings  may  be  set 
up.  For  non-planar  intersections,  limit  strings  have  three  elements, 
not  two.  For  basic  parameter  limits,  the  first  element  of  the  limit 
string  is  one  of  the  limits,  the  second  element  is  the  other  limit,  and 


the  third  element  is  the  same  as  the  first. 


9.2  Tracing  the  Edge  to  Find  Bounded  Parameter  Limits 


One  next  traces  the  edge  within  the  basic  parameter  limits.  The 
limit  string  is  <tj,  t2 , t3>,  with  tj  = t3.  First,  one  traces  from  tj 
to  t2,  and  then  one  traces  in  the  reverse  direction  from  t2  to  t3  = t3 . 
The  two  traces  use  opposite  "radical  signs"  in  the  solution  of  eq.  ( 37 ) - 

In  any  trace  from  t^  to  t^  , if  the  parameter  is  increasing 
< ti+1K  one  solves  eq.  (37)  with  a positive  radical  sign 

s = (-b+Zb^-Uac  )/2a  (kl) 

If  the  parameter  is  decreasing  (t^  > t^+^),  one  solves  it  with 
a negative  radical  sign 

s = (-b-/tZ-kac ) /2a  (42) 

This  implicitly  gives  the  radical  sign,  without  it  having  to  be 
specified  explicitly.  The  radical  sign  is  the  "Boolean  parameter" 
mentioned  earlier. 

As  one  is  tracing  the  curve,  one  looks  for  parameter  bounds  just 
as  in  section  7.3.  If  the  point  corresponding  to  t’  = t2  is  in  bounds, 
there  is  a bounded  parameter  limit  string  with  t2  as  its  second  element. 
All  other  bounded  parameter  limit  strings  would  have  two  limits.  Since 
limit  strings  for  non-planar  surfaces  should  have  three  elements,  the 
second  element  may  be  repeated.  Whenever  t^  = t^+^,  there  is  no  trace 
performed. 
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X.  LIMBS 


10.1  Defining  Limbs 

Since  quadric  surfaces  are  generally  not  flat,  it  is  likely  that, 
from  a particular  viewpoint,  the  quadric  surface  may  "fold  in  back  of 
itself".  The  locus  of  points  where  this  happens  is  the  limb.  This 
terminology  is  due  to  [Comba  1968].  A virtual  edge  is  a segment  of  a 
limb. 

More  precisely,  the  limb  is  the  locus  of  points  on  a surface  where 
the  normal  to  the  surface  is  perpendicular  to  the  line-of-sight . These 
points  must  satisfy  both  these  conditions: 


q(x,y,z)  = 0 and 


p = s.grad  q = 0 


where  s is  the  line-of-sight  vector  from  the  object  point  to  the  view- 
point; and  grad  q is  the  surface  normal,  which  may  be  expressed  as  a 
column  vector: 


I 2qjX  + q4y  + qfiz  + q? 
grad  q = 2q2y  + q4x  + qsz  + q8 


2q3z  + q6x  + q5y  + q? 


10.2  Orthographic  Projections 

For  simplicity,  let  us  assume  that  the  object  has  already  been 
translated  and  rotated  so  that  it  is  now  in  the  picture  space.  We  shall 
call  the  axes  x,  y,  and  z.  For  othographic  projections,  the  viewpoint  is 
considered  to  be  on  the  x-axis  at  infinity.  The  line-of-sight  vector  is, 
therefore,  along  the  x-axis  and  its  normalized  value  is  (100).  From 
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eq.  (1*3)  we  have 


p = s.grad  q =2qx+qy+qz+q 

14  6 7 


= 0 


(1*5) 


This  is  a plane.  It  is  called  the  "polar  plane  of  the  surface".  It 
does  not  exist  if  q = q = q = 0.  (A  paraboloid  whose  axis  is  the 

1 4 6 

x-axis  is  an  example  of  a surface  without  a polar  plane.) 


10.3  Perspective  Projections 

Woon  did  not  discuss  polar  planes  for  perspective  drawings,  and  it 
happens  that  they  do  require  a different  equation  than  orthographic 
drawings  do.  Assume  that  the  viewpoint  is  on  the  x-axis,  a distance  D 
from  the  origin.  We  make  the  simple  transformation  x'  = x - D,  putting 

the  viewpoint  at  the  x'-y-z  origin.  The  line-of-sight  vector  thus 

. . 

becomes  s = (-x'  -y  -z). 

Therefore,  the  equation  for  the  polar  surface  is 


p = s.grad  q 

= -2q'  x2-2q  y2-2q  z2-2q'  xy-2q  yz-2q'  zx'-q'  x'-q  y-q  z=0  (1+6) 

123  45  6 789 

By  substituting  twice  the  value  of  q(x,y,z)  from  eq.  (l),  we  have 


p = -2q’(x',y,z)  + q'  x ' + q y + q z + 2q  =0 

7 8 9 0 


(1*7) 


Since  all  points  on  the  limb  must  be  on  the  surface,  with 
q'(x',y,z)  = 0,  the  above  equation  reduces  to: 


p = q'  x + q y + q z + 2q  , 

v 7 8 9 0 


(1*8) 


where  p^  = 0 is  the  equation  of  the  "virtual  polar  plane". 
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10.  U Point  Orientation 

The  function  p(x,y,z)  has  another  use  besides  acting  as  the  equa- 
tion of  the  polar  plane.  It  also  indicates  the  orientation  of  a point 

*► 

on  the  surface.  If  p > 0,  the  angle  between  s and  the  normal  is  less 
than  90°,  and  so  the  point  is  front-oriented  and  potentially  visible. 
If  p < 0,  then  the  angle  is  greater  than  90°  and  the  point  is  back- 
oriented  and,  therefore,  invisible.  Of  course,  if  p = 0,  then  the 
angle  is  90°,  and  the  point  is  "orthogonally"  oriented  and  is  on  a 
limb. 

For  perspective  drawings  the  function  p^  may  be  used  in  place  of 
p,  since  all  points  that  ever  need  to  be  tested  by  the  value  of  p are 


on  the  surface  anyway. 
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XI.  HIDDEN-LINE  DETERMINATION  AND  PICTURE  DRAWING  | 

This  algorithm  is  primarily  desgined- to  handle  QSICs,  with  hidden- 
line  elimination  a secondary  consideration.  The  methods  described  here 
may  be  considered  "brute-force''  methods,  but  they  should  work  well. 

11.1  Preliminary  Processing 

Each  quadric  patch  is  considered,  and,  if  it  contains  a limb  within 
its  bounds,  it  is  divided  into  two  or  more  faces  at  the  limb.  If  there 
is  no  limb,  then  the  patch  is  a single  face. 

Each  face  is  tested  to  see  whether  it  is  a front-face  or  a back- 
face.  As  Woon  pointed  out,  [Woon  1970,  p.  26],  a face  has  the  same 
orientation  as  any  point  on  it. 

For  an  edge  to  be  visible  at  least  one  of  the  faces  adjacent  to  it 
must  be  a front-face.  This  algorithm  does  not  use  Loutrel's  complete 
classification  system  [Loutrel  1970]  because  it  is  dependent  on  whether  a 
surface  intersection  is  obtrusive  or  intrusive,  and  a QSIC  may  be  both. 

The  parameter  limit  strings  of  an  edge  may  have  to  be  changed.  If 
an  edge  is  entirely  the  intersection  of  two  back-faces,  it  will  not  be 
displayed  at  all  in  a hidden-line  eliminated  drawing.  If  it  is  the 
intersection  of  two  quadric  patches  with  limbs,  some  of  it  may  be  the 
intersection  of  two  back-faces  and  be  hidden,  while  the  rest  of  it  may 
be  potentially  visible.  The  limits  must  be  appropriately  reset  so  that 
only  potentially  visible  parts  are  included. 

At  this  point,  two  comments  are  in  order.  First,  limbs  will  now 
be  considered  to  be  like  any  other  intersection,  with  parameter  limits 
and  coefficients.  Second,  the  following  discussion  is  for  orthographic 


A 
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drawings  only.  For  perspective  drawings,  one  must  substitute  "y/x" 
for  "y"  and  "z/x"  for  "z",  wherever  these  appear. 

For  each  potentially  visible  edge,  the  y-  and  z-extrema  should 
be  computed.  (The  x-extremum  may  also  be  computed.)  This  may  be  done 
by  tracing  through  an  edge,  looking  for  places  where  each  coordinate 
has  an  extremum.  This  will  happen  when  the  sign  of  the  increment  of 
the  coordinate  changes.  This  method  is  like  looking  for  "dips". 

For  each  front-face,  the  extrema  of  the  extrema  of  the  associated 
edges  are  taken.  Thus  we  have  both  edge  extrema,  and  face  extrema. 

11 . 2 Tracing  an  Edge 

First,  for  each  edge  one  must  find  those  front  faces  which  may 
hide  all  or  part  of  the  edge.  This  is  done  by  using  the  edge  extrema 
and  the  face  extrema.  A face  may  hide  an  edge  only  if  the  y-extrema 
and  the  z-extrema  overlap.  (If  x-extrema  are  calculated,  the  nearer 
x-extrenum  of  the  face  must  be  closer  to  the  viewpoint  than  the  further 
x-extremum  of  the  edge  for  that  edge  to  be  potentially  hidden  by  the 
face.  ) 

One  is  now  ready  to  trace  the  edge,  simultaneously  detecting  which 
parts  are  hidden,  and  drawing  the  edge.  As  soon  as  one  is  sure  that  the 
vector  from  one  point  to  the  next  is  completely  unhidden  (hidden),  one 
may  draw  (skip)  that  vector.  Thus  both  hidden-line  detection  and  edge 
drawing  occur  as  parts  of  the  same  process. 

There  are  two  methods  of  scanning  an  edge  and  they  are  outlined 
below.  Either  may  be  used. 
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The  "tight"  scan  uses  small  vectors.  As  each  point  is  scanned, 
one  computes  the  line  of  sight  from  that  point  (test  point)  to  the 
vantage  point.  One  then  computes  the  intersection  of  this  line  with  the 
quadric  surface  associated  with  each  front  face  which  may  hide  the  point. 
This  is  a second-order  equation,  similar  to  eq.  (37).  If  there  are  no 
roots,  the  point  is  not  hidden  by  that  face.  If  there  are  two  roots, 
one  of  them  corresponds  to  a back-face,  and,  since  a back-face  cannot 
hide  a point,  this  root  is  discarded. 

One  takes  the  (remaining)  root,  and  checks  to  see  if  the  correspond- 
ing point  is  within  bounds  for  that  surface  patch.  If  it  is,  the  point 
is  hidden  if  the  intersection  point  is  between  the  test  point  and  the 
view  point.  If  the  intersection  point  is  out  of  bounds,  or  is  further 
from  the  view  point  than  the  test  point,  then  the  test  point  is  tested 
against  the  next  face.  If  no  face  hides  the  point,  then  it  is  visible. 

If  one  point  is  hidden  and  one  of  its  neighbors  is  not,  then  there 
is  some  point  between  them  where  the  edge  becomes  hidden.  By  using 
interval-halving,  one  may  locate  this  point. 

This  method  uses  many  small  vectors  to  insure  that  a small  hidden 
(unhidden)  section  of  the  edge  is  not  considered  visible  (hidden).  The 
disadvantage  of  this  method  is  that  many  points  have  to  be  inspected  for 
each  edge. 

The  "loose"  scan  uses  larger  vectors  but  needs  to  store  a good  deal 
of  information  about  each  point.  This  is  all  information  that  would 


have  to  be  computed  anyway. 
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For  each  front  face  which  could  hide  part  of  the  edge,  a list  called 
the  "state  list"  is  made  of  the  values  of  the  following  functions 
t l)  The  quadratic  discriminant  {V  = b2  - Uac)  of  the  second-degree 

equation  of  the  face's  surface  along  the  line  of  sight,  and 
2)  The  bounding  surface  function  q,  (x.,y.,z.)  of  each  bounding  surface 

Dill 

of  the  natch  on  which  the  face  is  located,  where  x. , y. , and  z.  are  the 

11  i 

coordinates  of  the  intersection  point. 

If  there  is  a change  of  sign  of  any  function  in  any  state  list,  then 
it  is  possible  that,  at  some  intermediate  point,  the  edge  becomes  hidden. 
One  may  use  interval  halving  the  locate  this  point.  It  is  also  helpful 
to  look  for  "dips"  in  these  functions. 

The  "loose"  scan  gives  an  effective  resolution  which  is  much  better 
than  the  vector  length  actually  used. 

11.3  Vector  Length 

In  using  the  "loose"  scan,  one  may  use  a variable  vector  length. 

If  p and  q are  two  adjacent  vectors,  then  |p|x|q|  = |p||q|sin  pq.  By 
trying  to  equalize  |p|x|q|  one  should  get  an  optimal  balance  between 
having  fewest  vectors  (requiring  least  processing  and/or  storage  space) 
and  making  sure  no  angles  are  too  sharp.  Using  this  equalization,  the 
parts  of  a curve  with  the  greatest  curvature  will  have  the  shortest 
vectors  and  the  largest  inter-vector  angles,  while  nearly  straight  curves 
would  have  long  vectors  and  small  angles.  However,  because  of  the 
requirements  of  the  algorithm  in  looking  for  sign  changes  and  "dips", 
each  edge  should  have  at  least  two  or  three  vectors.  It  might  also  be 
good  to  have  an  upper  limit  on  vector  length. 

---  - ■ _ 


r 
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XII . CONCLUSION 

The  chief  advantages  of  this  algorithm  are  (l)  it  is  a complete 
quadric  surface  algorithm,  allowing  use  of  all  real  quadric  surfaces; 

(2)  it  allows  "smooth"  intersections  to  approximate  higher-order  surfaces, 
and  (3)  it  may  be  fast  enough  to  be  used  in  conjunction  with  a shading 
algorithm  such  as  outlined  in  [Phong  1975 ] - 

Since  this  algorithm  distinguishes  limbs,  sharp  edges,  and  smooth 
edges,  each  can  be  appropriately  handled  by  the  shading  algorithm.  This 
can  include  explicit  handling  of  the  mach  band  effect,  specular  reflec- 
tion, and  transparency  [Metelli  197^]. 
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APPENDIX  1 - PROOF  OF  THEOREM 

In  this  report  it  has  been  assumed  that  all  QSICs  lie  either  in 
singular  quadric  surfaces  or  else  in  hyperbolic  paraboloids.  The  follow- 
ing is  a proof  of  this  conjecture. 


DEFINITIONS 

A quadric  surface  is  positive , singular , or  negat ive  if  the  deter- 
minant of  its  discriminant  matrix  is  positive,  zero,  or  negative  respec- 
tively. A pencil  of  quadric  surfaces  is  singular  iff  it  contains  at 
least  one  singular  quadric  surface,  and  is  non-singular  otherwise.  A 
surface  intersection  (QSIC)  is  singular  if  the  pencil  of  the  two  inter- 
section surfaces  is  singular,  and  is  non-singular  if  the  pencil  is 
non- singular. 

A congruence  transformation  is  a non-singular  transformation  such 
that  the  upper-left  3x3  submatrix  of  the  transformation  matrix  (called 
the  rotational  part)  is  orthogonal,  and  the  last  column  of  the  matrix  is 
(000  i).T 

A scaling  transformation  is  a non-singular  transformation  of  the 

3 

form: 


I a 0 0 0 
0 6 0 0 
0 0 Y 0 
0 0 0 6 


with  a,  3,  Y>  <$  f 0. 


A valid  transformation  is  the  product  of  congruence  and  scaling 
transformations . 

Notice  that  valid  transformations  do  not  change  the  type  of  surface 
the  discriminant  describes,  although  scaling  transformations  may  distort 


such  a surface. 


DEFINITION: 

A PARA  is  a quadric  surface  with  a singular  sub-discriminant. 

PARAs  include  all  singular  quadric  surfaces  except  cones,  as  well  as 
elliptic  and  hyperbolic  paraboloids. 

LEMMA  1: 

The  pencil  of  any  two  quadric  surfaces  contains  a PARA. 

Proof:  Suppose  we  have  two  quadric  surfaces,  P and  Q. 

[Case  1:  Either  P or  Q is  a PARA.  The  hypothesis  is  satisfied. 

[Case  2:  Neither  P nor  Q is  a PARA.  The  equation: 

det(Ru(a))  = detCQ^-aP^)  = 0 is  of  the  form: 

-det ( P )a3  + K a2  - K a + det(Q  ) = 0, 
u 2 l u 

with  K and  K being  the  sums  of  the  determinants  of  combinations  of  the 
2 l 

columns  of  P and  Q . Because  det(P  ) # 0 and  det(Q  ) f 0,  this  equation 
u u u u 

is  of  third  order  and  no  less.  Since  a third  order  equation  must  have 
at  least  one  real  root,  there  is  at  least  one  value  of  a such  that 
det(R(a))  = 0.  This  R(a)  is,  by  definition,  a PARA. 

LEMMA  2: 

Any  pencil  of  a positive  quadric  surface  and  a negative  quadric 
surface  is  singular. 

Proof:  The  equation  det(R(a))  = det(Q-aP)  = 0 is  equivalent  to 

det(P)a4  -Ka3+Ka2-Ka+  det(Q)  = 0,  where  K , K , and  K are 

3 2 1 3 2 1 

sums  of  determinants  of  combinations  of  columns  of  P and  Q.  If  P is 

positive  and  Q is  negative,  then  det(R(0))  = det(Q)  < 0;  while  for 

a -►  ± °°,  det(R(a))  -*•  + 00  > 0.  Because  R(a)  is  a continuous  function, 

J 


there  must  be  at  least  two  values  of  a,  one  positive  and  one  negative, 
for  which  R(a)  is  singular. 


-A-3- 


LEMMA  3: 

One  may  apply  any  valid  transformation  to  any  pair  of  quadric  sur- 
faces without  affecting  the  roots  of  the  equations:  det(Q-aP)  = 0 and 

det(Q  -BP  ) = 0. 
u u 

Proof:  Suppose  we  have  the  valid  transformation  S and  its  rotational 

part,  S . Both  S and  S are  non-singular.  If  Z = S \ then  Z = S \ 
u u * u u 

and  both  S and  Z are  non-singular, 
u u 

Det(SQE-aSPZ)  = det(S(Q-aP)I)  = det(S)  det(Q-aP)  det(E) 

= det(S)det(£)det(Q-aP)  = det(Q-aP).  Similarly  for  and-  8. 

LEMMA  k: 

The  arbitrary  elliptic  paraboloid  P and  the  arbitrary  quadric  sur- 
face Q may  be  expressed,  in  some  transformed  space,  as: 


Q = 


Proof:  First,  let  us  select  the  axes  such  that  the  z-axis  is  the  main 

axis  of  the  paraboloid  P,  the  origin  is  located  on  P,  and  the  x-  and 
y-axes  are  the  axes  of  the  ellipse  which  the  cross-section  of  P.  This 
puts  the  surfaces'  discriminants  in  the  form: 


lA 

D 

0 

\ 

G ) 

0 

0 

x\ 

D 

B 

E 

H 

0 

l 

0 

Y 

22 

0 

E 

C 

J i 

P = 

0 

0 

0 

1 

H 

J 

K ) 

U 

Y 

1 

0/ 

Q"  = 


P”  = 


/ S 0 0 X \ 
0 T 0 Y 


0 0 0 Z 


\ 


X Y Z 


0 / 


r 


1 


If  S is  negative,  reverse  the  sign  of  every  element  of  P". 

Now,  use  the  following  scaling  transformation  on  both  P"  and  Q": 


This  gives  us: 


/ A 

D 

F 

A 

0 

0 

X \ 

D 

B 

E 

H 

| P’  = 

f o 

1 

0 

Y 

F 

E 

C 

J 

I 

1° 

0 

0 

1 

\G 

H 

J 

* ) 

\ X 

y 

l 

0 / 

We  now  apply  the  following  rotation  (which  is  a congruence  trans- 
formation) to  both  P'  and  Qf: 

L -M  0 
M L 0 
0 0 1 
\0  0 0 


Using  this  transformation  yields: 


j A D 0 G \ 

/ 1 0 0 X ' 

D 3 E H 

0 E C J 

lG  H J K j 

P = 

0 1 0 Y 

0 0 0 1 

ix  Y 1 0 ) 

0 \ 
0 
0 
1 


/ 


with  L2  + M2  = 1 

. L _ E 

and  M " F • 


/ 
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LEMMA  5: 

The  pencil  of  two  elliptic  paraboloids  is  always  singular. 

Proof:  Suppose  we  have  two  elliptic  paraboloids  P and  Q,  in  the  form 

given  in  LEMMA  There  are  some  conditions  that  are  imposed  on  the 
elements  of  Q,  as  follows: 


h 

0 

0 

/A 

D 

0 

G 

0 

1 

0 

Y 

Q = 

' D 

B 

E 

H 

P = 

0 

0 

0 

1 

0 

E 

C 

J 

Y 

1 

0 ) 

f 

iG 

H 

J 

K / 

det(P  ) = 0 
u 

T (P)  > 0 
2 

The  two  conditions  for  Q are 
2 


such  that:  det(Q^)  = 0, 

T (Q)  > 0. 

2 


2 _ 

AB  + BC  + AC  - D2  - E2  > 0. 


det(Qu)  = ABC  - CD2  - AE2  = 0; 


T (Q) 
2 


[Case  1:  C= 0:  det(Q^)  = 0 only  if  A = 0 or  E = 0. 


[Case  la:  C = E = 0 

Q = 

Obviously,  det  (Q-JP)  = 0. 
[Case  lb:  A = C = 0. 


/ A D 0 G 
D B 0 H 


0 0 0 J 

G H J K 


/0  D 0 G 
D B E H 
0 E 0 J 
G H J K 


Now,  T (Q)  = - D2  - E2  £ 0.  However,  one  of  the  conditions  for  Q to  be 
2 

an  elliptic  paraboloid  is  that  T >0.  Thus,  case  lb  never  occurs. 

2 


How,  det(Q  ) = A3  - AE2  - D2  = 0,  so  D = ±/A ( 3-Sz ) . 
u 

For  D no  be  real,  we  must  nave: 

A(3-E2)  2l  0 (condition  a) 

We  also  have:  T (Q)  = A3  + A + B - D2  - E2  > 0. 

2 

Since:  D2  = AB  - AS2 , ve  have: 

T (Q)  = A + 3 - E2  + A.S2  (condition  b) 

2 

new,  taking  the  pencil  of  P and  Q:  R(a)  = (Q-aP);  solving 
det(R  (a))  = det(<3  -a?  ) = 0,  we  get: 

U 'll 

- det{?  )a3  + a2  - (A+3-E2 )a  + det(Q  ) = C. 
u u 

Since  det(?  ) = det(Q  ) = C,  when  we  reject  the  root  a = 0,  we  have 
''-l  ^ 

the  remaining  root:  a = A + 3 - £2 . For  this  a,  F.(a)  is  a PARA,  by 

■"i  ^ f*  "L  n i **  i.  o ^ 

[Case  2a:  If  R(a)  is  a singular  PARA,  then  the  hypothesis  is  satisfied. 

[Case  2b:  If  R(a)  is  not  a singular  PARA,  it  must  be  a paraboloid, 

either  elliptic  or  hyperbolic.  We  will  assume  that  it  is  elliptic. 

Then  I (S;)  > 0. 

2 

Substituting  the  root  of  a as  above,  we  have: 

I E2-3  D 0 \ 

R (a)  = (Q  -a?  ) = I D E2-A  E 1 

U -1  u I 

0 E 1 . 


r 


■ 


Ta> 

ir.g  the  condition  that  T (R(a))  > 0: 

2 

(Taking  D2 

= AS  - AE2! 

) 

T = (E2-A)(E2-3)  +2E2  - (A+3) 
2 

- A3  + 

AE2 

> 0; 

T = E2(E2-3)  + E2  - A - 3 > C 
2 

i;  or 

-T  = E2(5-E2)  + A + B - E2  < 
2 

0. 

Row,  taking  S = 3 - E2 , the  three 

conditions 

are : 

A0  > 0 

(condition 

a) 

A + 3 + AE2  = A(l+E2)  + 3 > 0 

( condition 

b) 

E23  + A + 3 = A + 3(1+E2  ) < 0 

(condition 

c) 

By 

(a),  if  3 < 0,  then  A < 0;  if  A < 

0,  then 

3 < 

C; 

if  3 > 0,  then  A > 0;  if  A > 

0,  then 

3 > 

0. 

(b),  if  3 £ 0,  then  A > 0;  if  A £ 

0,  then 

3 > 

0. 

3y 

!o),  if  8 > 0,  then  A < 0;  if  A > 

0,  then 

3 < 
. / \ 

0. 

These  conditions  are  inconsistent;  hence  R(a)  cannot  he  an 
elliptic  paraboloid  and  roust  he  hyperbolic  paraboloid.  Because  this  is 
a pencil  of  tvo  negative  elliptic  paraboloids  and  a positive  hyperbolic 
paraboloid,  by  LEMMA  2,  the  pencil  and  its  associated  Q3IC  are 
singular . 

LEMMA  t : 

The  pencil  of  tvo  negative  quadric  surfaces  is  singular. 


Otherwise,  the  PARA  is  an  elliptic  paraboloid.  Let  us  call  the 


PAPA  ?,  and  choose  one  of  the  two  negative  quadric  surfaces  Q.  If  0, 
another  elliptic  paraboloid,  then  the  pencil  is  singular  by  LEMMA  5. 
Otherwise,  Q must  be  an  ellipsoid  or  a hyperboloid  of  two  sheets,  wi 
det(Q  ) 7 0.  l 

3 y LEMMA  L,  the  discriminants  of  P and  0 may  be  expressed: 


/ A D 0 G \ 

1 0 0 X \ 

D 3 3 H 

1 ? = 

0 1 0 Y 

0 S C J 

1 

0 0 0 1 

\G  H >J  K / 

l X Y 1 0 J 

The  solution  of  the  equation:  det(R  (a))  = det(Q  -a?  ) = 0 is: 

u u u 

-det (?_  )a2  + Ccr  - (3C-E2+AC)a  + det(Q  ) = 0;  or 

U U. 

Ca2  - [ (A+3)C-E2 ]a  + [ (A3-D2 )C-AE2 ] = 0 
[Case  1:  C = C.  Equation  (Al)  gives  a = A. 


/ 0 D 0 , 


?.  (A)  * U -A?  ) » D B-A  E 

u u u 


T ( R ( A ) ) = -D2-!2  < 0. 
2 — 


If  D = a = 0,  then  Rank  (R  (A))  is  as  most  one.  Since  the  rar.l-; 

u 

the  discriminant  cannot  exceed  the  rank  cf  the  sub-dsicriminant  by  r 
than  two.  Rank  (R(A))  <_  3,  so  R(a)  is  singular,  and  so  is  the  pencil 
the  QSIC. 

If  D ? 0 or  E ? 0,  then  T (^.(a) ) < 0,  so  R^A})  is  a hyperbolic 

2 

paraboloid,  and,  by  LEMMA  2,  the  QSIC  is  singular. 

[Case  2:  C f C.  If  we  normalize  so  that  C = 1,  equation  (Al)  beccir 

a2  - [(A+3)-E2]a  + [A3-D2-AE2]  = 0.  Taking  the  quadratic  form: 


-A-9- 


acr  + ba  + c = 0,  we  get  the  quadratic  discriminant:  = b2-hac. 

For  this  equation,  V = (A-B+E2)2  + LD2  > 0. 


[Case  2a:  If  = C,  then  D = 0 and  S2  = 


A.  There  is  cne  double 


root,  a = A,  and  Rank  (R^(A) ) <_  1,  so  Rank  (R(A))  <_  3 , so  R{a)  is 
singular,  and  so  is  the  pencil. 

[Case  2b:  If*?  > 0,  then  there  are  two  roots  in  a.  Both  R (a)  are 
PARAs.  If  one  is  singular,  the  hypothesis  is  satisfied.  If  one  is  a 
positive  hyperbolic  paraboloid,  then  the  pencil  is  singular  by  LEMMA  2. 
If  cne  is  an  elliptic  paraboloid,  the  pencil  is  singular  by  LEMMA  5. 
THEOREM : 

The  intersection  of  two  quadric  surfaces  either  is  singular,  lies 
in  a hyperbolic  paraboloid,  or  both. 

Proof : If  cne  of  the  surfaces  is  singular,  the  hypothesis  is  satisfied. 

If  one  is  positive  and  cne  negative,  the  pencil,  and,  therefore,  the 
intersection,  is  singular  by  LEMMA  2.  If  both  are  negative,  then  the 
pencil  and  QSIC  are  singular  by  LEMMA  6.  If  both  are  positive,  then  by 
LEMMA  1,  the  intersection  is  contained  in  a PARA.  If  the  PARA  is 
negative  or  singular,  the  intersection  is  singular  by  LEMMA  2,  or 
directly.  Otherwise,  the  PARA  must,  be  a hyperbolic  paraboloid. 
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